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Abstract 

We study well-posedness of sweeping processes with stochastic perturbations gener¬ 
ated by a fractional Brownian motion and convergence of associated numerical schemes. 

To this end, we first prove new existence, uniqueness and approximation results for de¬ 
terministic sweeping processes with bounded p-variation and next we apply them to the 
stochastic case. 
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1 Introduction 

In the present paper we study well-posedness of some variants of the so-called sweeping 
process introduced by Moreau in the early 70s with motivation in plasticity theory. In his 
original formulation the sweeping process coincides with a first order differential inclusion of 
the form 

r f(t)€iV(a;x(t)), 

I x(0) = xo € Co, (1.1) 

[ x{t) G Ct, 

where Ct is a given convex moving set and N(Ct',x{t)) is the inward normal cone to Ct at 
point x{t) (see [351 [HHl EZ] )■ Many attempts have been made to generalize Moreau’s results 
to larger class of moving sets or more general than dEU) differential inclusions containing 
deterministic or stochastic perturbations. For instance, sweeping by prox-regular moving 
sets instead of convex sets was considered by Colombo and Goncharov [12], Benabdellah 
[H, Thibault [M], Colombo and Monteiro Marques m- The study of sweeping processes 
with perturbations was introduced by Castaing, Due, Ha and Valadier [9] and Castaing amd 
Monteiro Marques m- The interest in the theory of sweeping processes comes from the fact 
that it has numerous practical applications in nonsmooth mechanics, analysis of hysteresis 
phenomena, mathematical economics and in the modeling of switched electrical circuits (see, 
e.g., the monographs by Acary, Bonnefon and Brogliato [1], Drabek, Krejci and Takac |16] . 
Monteiro Marques |34j and the references therein). 
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In our paper we study sweeping processes with stochastic perturbations. This problem 
was considered earlier by Colombo SIB] and recently by Bernicot and Venel [S]. In the 
last paper the authors give conditions ensuring well-posedness of d-dimensional stochastic 
differential inclusions of the form 

dXt G fit, Xt)dt + g{t, Xt)dBt + NiCt-, Xt), 

Xq = xq € Co, (1.2) 

G Ct, 

where Ct is a given prox-regular moving set and B = {Bt}t£M+ is a standard Brownian 
motion. To do this, in proofs they combine the methods of deterministic sweeping process 
theory with the methods of stochastic differential equations (SDEs) with reflecting boundary 
conditions. The use of the methods of SDEs is possible, because one can observe that (II.2p 
is equivalent to the SDE with reflecting boundary condition of the form 

Xt = xo+ f f{s,Xs)ds + 

Jo 

where the integral with respect to B is the classical stochastic integral. 

By a solution to Oi we mean a pair {X,K) consisting of a process X = {X}^g]K+ 
such that Xf G Ct and the process K = {Kt}t£M.+ , called regulator term, such that dKt G 
N{Ct;Xt) in appropriately defined sense. Equation (II.3h was firstly investigated by Sko- 
rokhod [H] for Ct = [0,oo), t G M^. Extensions of Skorokhod’s results to larger class of do¬ 
mains was studied for instance by Tanaka [50], Lions and Sznitman [32], Saisho [33], Dupuis 
and Ishi |20j . Slomihski |46j and Rozkosz [40] . Equations of the form (jl.3p also have many 
applications, for instance in queueing systems, seismic reliability analysis and finance (see, 
e.g., [3[EU[28[HI] and the references therein). Solutions of (II.3p are often called solutions of 
Skorokhod’s SDEs or of the Skorokhod problem. 

In our paper a stochastic perturbation is generated not by a standard Brownian motion 
but by a fractional Brownian motion (fBm) B^ = with Hurst index H > 1/2, i.e. 

by a continuous centered Gaussian process with covariance 

= i(tP + - \t2 - tip^), h,t2 G M+. 

It is well known that B^ is not a semimartingale and therefore the classical stochastic inte¬ 
gration theory for semimartingales cannot be applied. However, B^ has A-Holder continuous 
paths for all A G (0, Lf), which allows one to define the pathwise Riemann-Stieltjes integral 
with respectto fBm (see, e.g., [isliiasi])- The theory of SDEs without reflecting boundary 
condition driven by B^ with the pathwise Riemann-Stieltjes integral is at present quite well- 
developed. General results on existence and uniqueness of solutions one can find in Nualart 
and Ra§canu [39]. The viability property for such equations is considered in details in Giotir 
and Ra§canu [H]- 

Let / : M'’* —>■ g : ^ he measurable functions and B^ be a d-dimensional 

fBm. Our main purpose is to study d-dimensional SDE with reflecting boundary condition 
of the form 

Xt = Xo+ f fis,Xs)ds+ f gis,Xs)dB^+ Kt, t G M+, (1.4) 

Jo Jo 


gis, Xg) dBg + Kt, t€ 


(1.3) 
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where the integral with respect to is the pathwise Riemann-Stieltjes integral and Ct = 
[Lt, Ut] = [//] C is a moving convex set (Here L\ < Ul, t ^ M+, i = 1,... ,d). 

We also study some generalizations of (|1.4n . Clearly, (HI is equivalent to sweeping process 
of the form (jl.2p with stochastic perturbation generated by fBm. The integral form (jl.4p is 
however more convenient because in general the process K need not be of bounded variation 
and therefore the use of the differential dKt would require additional explanations. 

In the recent paper by Ferrante and Rovira [2Tj the special case of HID with Ct = [0, ooY 
was considered. Using quite natural in the context of SDEs driven by methods based on 
A-Holder norms they gave conditions ensuring the existence of solutions and their uniqueness 
for some small time interval. Some global uniqueness results for (jl.4p with time homogenous 
coefficients f,g and Ct = xf^^[Ll,oo) were proved in Falkowski and Slomihski [23], where 
in contrast to [23| the p-variation norm is used. In the present paper we also use techniques 
using the p-variation norm. It is worth noting that we do not assume the so-called “interior 
ball condition”, which in our case means that there is r > 0 such that Ul — L\ > r, t € M"*". 
We even allow that LI = UJ:. Unfortunately, we are not able to extend our methods to more 
general moving convex sets or prox-regular moving sets (we think that it is not possible apart 
from the case of functions g depending purely on time). 

As a matter of fact, in the present paper we study more general than (11.411 equations in 
which the driving processes may have jumps, that is equations of the form 

f {s, Xg-) dAg + f g{s, Xg-) dZg + Kt, t € M'*'. (1-5) 

Jo 

where A is a one-dimensional cadlag process with locally bounded variation and Z is a, d- 
dimensional cadlag process with locally bounded p-variation for some 1 < p < 2 (note that 
B^ has locally bounded p-variation only for p G (1/if, oo)). 

The paper is organized as follows. 

In Section 2 we consider the deterministic extended Skorokhod problem x = y + k as¬ 
sociated with a cadlag d-dimensional function y (i.e. y G D(M'’',M'^)) and time dependent 
barriers l,u € D(R+,M'^) such that I < u, which means that It < ut, t ^ lo < yo < uq. 

We show that for fixed I, u the mapping y i—>■ (x, k) is Lipschitz continuous in the p-variation 
norm. It is worth noting here that in m Remark 3.6] it is observed that y i->- (x, k) is not 
Lipschitz continuous in the A-Holder norm and that for that reason in [23| the authors were 
not able to obtain global uniqueness. 

In Section 3 we consider a deterministic version of (11.511 . We give conditions ensuring the 
existence and uniqueness of solutions. In the proof we use an analogue of the Picard iteration 
method (we work in spaces equipped with the p-variation norm). Our assumptions on the 
coefficients f,g are similar to those considered in |39|. Since our integrators are cadlag with 
bounded p-variation and need not be A-Holder continuous, our theorem generalizes results 
from m even in the trivial case where Ct = R'^, t G R^. 

Section 4 is devoted to the approximation of deterministic solutions considered in Section 
3. We consider two methods of approximations. The first one is an easy to implement 
discrete-time method constructed in analogy with the classical Euler scheme (it is an analogue 
of the so-called “catching-up” algorithm introduced by Moreau to prove the existence of a 
solution to HID)- We prove the convergence of the scheme in the Skorokhod topology Ji 
(in the case of continuous data we obtain uniform convergence on compact subsets of R+). 
The second method uses stability of solutions of deterministic versions of (jl.5p with respect 
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to convergence of its coefficients. More precisely, we consider family of solutions with the 
coefficients fe,ge, e > 0 instead of /, g and such that —> /, g^ —> g' as e —)■ 0, which means 


K 


K. 


that fe,ge tend to f,g uniformly on compact subsets of We show that under some mild 
additional assumptions on fe,ge the associated solutions converge in the p-variation norm to 
the solution of equation with coefficients /, g. 

In Section 5 we apply our deterministic results to show the the existence and uniqueness 
of solutions of SDEs of the form m- To illustrate how our results work in practice we 
consider fractional SDEs ()1.4[) and its simple generalizations. We give conditions ensuring 
the existence and uniqueness of their solutions and show how approximate them by a simply 
to implement numerical scheme. 

Section 6 contains the proof of Theorem 12.21 

In the sequel we will use the following notation. M"*" = [0,oo), is the space of d x d 
real matrices A, with the matrix norm ||j 4|| = sup{|j4u|;u € |u| = 1}, where | • | denotes 

the usual Euclidean norm in B{0,N) = {x € |x| < A^}, N G M+. is 

the space of cadlag mappings x : M'*' ^ i.e. mappings which are right continuous and 
admit left-hands limits equipped with the Skorokhod topology Ji. For x € t > 0, 

we write xt- = lim^tiXs, Axt = xt - xt- and Vp{x)[a,b] = sup.,^ 1]?=! < +oo, 

where the supremum is taken over all subdivisions vr = {a = to < • • • < tn = b} of [a, b]. 
^p{x)[a,b] = ypix)[a,b] + \xa\, where Vp{x\a,b] = {vp{x)[a,b]Y^^, is the usual p-variation norm. 
Moreover, for simplicity of notation we write Vp{x)T = Up(x)[o,T]) yp{x)T = bp(®)[o,r] and 
yp{x)T = yp{x)]p Ty If X € D(]R+,M'^) then in the definition of p-variation Vp we use the 
matrix norm || • || in place of the Euclidean norm. We write x < x', x,x' G D(R+,R'^) if 


Xt < x't, t G R+, i = 1,... ,d. 


2 Main estimates 

We begin with recalling the dehnition of the extended Skorokhod problem with time depen¬ 
dent reflecting barriers introduced in m- Let y,l,u G D(R+,R'^) be such that I < u and 
^0 < 2/0 < ^ 0 - We say that a pair (x, k) G B(R+,R^'^) is a solution of the extended Skorokhod 
problem associated with y and barriers l,u (and we write {x,k) = ESP{y,l,u)) if 

(i) xt = yt + kt € [lt,ut], t G R+, 

(ii) ko = 0, k = {k ^,..., k'^), where for every 0 <t < q and i = 1 ,... ,d, 

kq — kl> 0, if x\ < u\ for all s G (t, g], 

kq — kl < 0, if X* > II for all s G (t, q], 

and for every t G R"*", Akl > 0 if xj < uj and Akl < 0 if xj > /J. 

In O Theorem 2.6] it is proved that for any y,l,u € B(R+, R*^) such that I < u and /q < 2/0 < 
uq there exists a unique solution {x,k) = ESP{y,l,u). 

Remark 2.1 (a) It is observed in |48] that instead of (ii) the following system of conditions 
can be considered: for every 0 < t and i = 1 ,... ,d such that inf^gj^,j](u® — /*) > 0 the 
function E has bounded variation on [t, q] and 

f (x* — ID dk\ < 0 and f (x* — u*) dk\ < 0 (2.1) 

J [t,q] J [t,q] 
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(Here Wg dvg denotes the integral over the closed interval [t, g], which is equal to wt^vt + 
Wg dvg, where Wg dvg denotes the usual integral over the half open interval (t, q]). Simple 
calculations show that the definitions from [6] and [38] are equivalent. 

(b) By (|2.1h . if ul > l\ then {x\ — ll)Akl < 0 and {xl — ul)Akl < 0. Consequently, if 
Akl > 0 then xl = II and if Ak\ < 0 then x\ = u\, i = 1,..., d. Therefore, for every t € R+, 

xt = max(min((xt_ + Ayt),ut), k) and h = max(min(/ct_, - yt),lt - yt), 

which means that xt is the projection of xt- + Ayt on the interval [ut,lt] and kt is the 
projection of kt- on the interval [ut — yt,lt — yt]- 

(c) In the classical Skorokhod problem it is assumed that the function k has bounded 

variation on each bounded interval [t, q], or, equivalently, k = — k^~\ where 

are nondecreasing right continuous functions with ko = k^'^ = k^ ^ = 0 such that 
increases only on {t‘,xl = Pt} and increases only on {t;xl = rt*t}, i = 1, ... ,d. If 

{x,k) = ESP{y,l,u) and mit<T{ut — k) > €t > 0, T & then A: is a function of bounded 
variation and {x,k) is a solution of the classical Skorokhod problem (see, e.g., (6] Corollary 
2.4]). 

The Lipschitz continuity of the mapping (y,l,u) i—>■ {x,k) in the supremum norm is well 
known. Let {x,k) = ESP{y,l,u), {x',y') = ESP{k',1',u'). By [48l Theorem 2.1], 

sup \xt - x't\ < 2sup \yt - y't\ + supmax(]Zt - Q, \ut - u't\) 
t<T t<T t<T 

and 

sup \kt - k[\ < sup \yt - y(] + supmax(]Zi - /(], \ut - rt(]). 
t<T t<T t<T 

From this one can deduce the following stability result for solutions of the extended Skorokhod 
problem in the topology Ji. Assume that [x'^^kP] = ESP{y"',E,u"'), n G N, {x,k) = 
ESP{y,l,u) and {y^,E,u'^) —^ {y,l,u) in D(R+,R3'^). Then 

{x^, k'^,y^, r, u”) ^ (x, k, y, I, u) in D(R+, R^^^) (2.2) 

(see |6l Theorem 2.6] or |48l Theorem 2.8]). Below we show that in the case of fixed barriers 
l,u the Lipschitz continuity of the mapping y (x. A:) also holds in the p-variation norm. 
We first consider the case d = 1. 

Theorem 2.2 Assume that y^,y‘^,l,u € B(R+,R) are such that Iq < yo,yQ < uq and I < u. 
Let {x^,k^) = ESP{y^,l,u), j = 1,2. Then for any T € R+ 

Vp{k^ - k^)T < Vp{y^ - y'^)T- 

Since our proof involves some technical one-dimensional arguments not associated with the 
rest of the paper, we defer the proof of Theorem 12.21 to Section 6. 

Remark 2.3 (a) The case p = 1 was studied earlier in |48l Theorem 2.14] (see also |43jL 

(b) In Ferrante and Rovira |24l Remark 3.6] it is observed that property stated in Theorem 
12.21 does not hold in A-Holder norm. 

(c) |48l Example 2.15] shows that it is not possible to omit the assumption that I = I' 
and u = u'. 
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Corollary 2.4 Assume y,y',l,u € are such that Iq < yo,yQ < uq. Let {x,k) = 

ESP{y, I, u) and {x', k') = ESP{y’, l,u). Then for any T € M+, 

Vp{x - x')t <{d+ l)Vp{y - y')T and Vp{k - k')T < dVp{y - y')T- 

Proof. By Theorem 12.21 

d d 

Vp{k - k')T < d^P-^'^/P{Y^Vp{P - k'^)T)^/P < d^P-^^/P{Y^Vp{y^ - y'%)^/P 

i=l i=l 

< (imaxPp(y* - 7 /'*)t < dVp{y - y')T- 

i 

Since Vp{x — x')t < Vpiy — y')T + Ppik — k')T, the proof is complete. □ 

Corollary 2.5 Assume y,l,h,u € are such that Iq < yo < uo, I < h < u. Let 

{x, k) = ESP{y, I, u). Then for any T € 

%{x)t <{d + l)Vp{y)T + dVp{h)T and Vp{k)T < dVp{y)T + dVp{h)T. 

Proof. Note that (/i,0) = ESP{h,l,u). By Corollary 12.41 

Vp{k)T < dVpiy - h)T < dVp{y)T + dVp{h)T, 

i.e. the second inequality of the corollary is satished. From the second inequality we imme¬ 
diately get the hrst one. □ 

3 Deterministic equations with reflecting boundary condition 

Let a G z,l,u G be such that Pi(a)T, Vp{z)T < oo for T G and 

I < u. We also assume that there is /i G B(R''',M'^) such that I < h < u and Vp{h)T < oo 
for T G R^. This additional assumption is indispensable to ensure that (x, k) = ESP{y, I, u) 
have bounded p-variation for any bounded p-variation function y (it is automatically satished 
if — lt)>eT>0,T£ R+, because in this case A: is a function of bounded variation). 

We consider equations with reflecting time-dependent barriers of the form 

Xt = xo+ f{s,Xs-)das+ g{s,Xs-)dzs + kt, f G R+, (3.1) 

Jo Jo 

where / : R+ x R*^ —^ R*^, g : R+ x R*^ —> are given functions, the integral with respect 

to 2 ; is a Riemann-Stieltjes integral and Zq < xq < uq. We recall that if rc G B(R''',M'^), 
2 ; G B(R+,R'^) are such that Vq{w)T < + 00 , Vp{z)T < -|-oo, T G R"*", where l/p+ l/q > 
1, p, g > 1, then the Riemann-Stieltjes integral f^Wg-dzs is well dehned (see, e.g., [HI). 
Moreover, it is well known that for any a < b, 

Pp{ 'a!s—dZs)[afi] ^ Cp^qVq{w)^a,b)^pi^)[a,b]j (3-2) 

J a 

where Cp^q = C(p~^ + Q~^) and f denotes the Riemann zeta function, i.e. ({x) = 1/^^- 
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Definition 3.1 We say that a pair {x,k) G is a solution of (|3.ip if Vp{x)T < oo 

for T G and {x,k) = ESP{y,l,u), where 



We will need the following conditions on /, g. 


(F) (a) There exists L > 0 such that 

\f{t,x)\ < L{1 + |x|), X G t G 
(b) For every N G there exists > 0 such that 

- f{t,y)\ < Ln\x - y\, x,y € 5(0,iV), t G M+. 

(G) (a) There exist /3 G (1 — 1/p, 1] and > 0 such that 


\g{t, x) - g{s, y)\ < C^{\t - s|^ + |x - y\), x,y G t, s G 


(b) g is differentiable in x and for every N G M+ there exist G {p—1,1] and Cn > 0 


such that 






Similar sets of conditions were considered in papers on equations without reflecting bound¬ 
ary condition driven by functions (processes) with bounded p-variation (see [I?! [29l |30l |33l 

I391I1I])- 

Remark 3.2 Note that under (G)(a) for every T G M"*" there exists C^’'^ > 0 such that for 
every t G [0, T] and x G 


\g{t,x)\ < + |x|), 

and for q=p\/ (1//3) and every w G B(M+,M‘^), 

%{g{;w))t < + C%{w)t + |5(0,xo)| < + Vp{x)t). 


(3.3) 


(3.4) 


Moreover, C^’’^ = C^{T^ + 1) + |5(0,0)|. 

We will approximate solutions of ()3.ip by using an analogue of the Picard iteration 
method. Set {x^,k^) = ESP{xq,1 , u) and for any n G N set 


xo + /o f{s, x'/J) das + /o g{s, x/J) dzs, 
ESP{y^,l,u), 



(3.5) 


where the integral with respect to z is the Riemann-Stieltjes integral. Note that (j3.5p is well 
defined if (F)(a) and (G)(a) are satisfied. Indeed, by Corollary 12.51 


Vp{x^)T <{d + l)\xo\+dVp{h)T, rGM+ 


(3.6) 
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and for any n G N, 


VJx^h <{d+ l)t4(y”)r + dVJh)T 


< (d + 1) 


\xo\ + Vp{f f{s,x'^_^)das)T + Vp{f g{s,x'^_^)dzsyj 
Jo Jo 


+ dVpih)T. 


Moreover, 


Vpi [ f (s, das)T < sup \f{s,x'^_y\Vi{a)T < Lil + Vpix"^ ^)T)Vi(a)r 

Jo s<T 

and by (|3.2p and (|3.4p . for g = p V (1//3) we have 

Vp{ [ g{s, x'^zy dzsh < CpMg{; x”-1))tFp(^)t < + Vp{x^-yT)Vp{z)T■ 

Jo 

Hence, in particular, Vp{x^)T < oo for n G N, T G M"*". In fact, under (F)(a) and (G)(a) we 
have 

supt^(x"')'r < oo, T G (3.7) 

n 

To check this, fix T G and set Cq = {d+ l)|xo| + dVp{h)T, Ci = {d+ 1) niax(L, C'p^gC'^’'^). 
Observe that by the above estimates for any t <T we have Vp{x^)t < Cq and 

Vpix'^y < Co + Ci(l + Vp{x^-yt){Vi{a)t + Vp{z)t), n G N. 

If we set ti = inf{t; Ci{Vi{a)t + Vp(z)t) > 1/2} A T then 

Vp{x^)t,-<Co + ^ + ^Vp{x^-%-, nGN, 
which implies that sup„ < 2(Co + 1/2). Since 

|Ax/J < \f{ti,x'^y})Aaty + \g{ti,xy~_})Azty + max(|A/tJ, |), 

it is clear that 

supVp{x"')ti < oo. (3.8) 

n 

Set tk = inf{t > 4-i;Ci(I/i(a)[t^_^^t] + Vp{z)[t^^_^^t]) > 1/2} AT, k >2. Modifying slightly 
the proof of (13.81) on can show that sup^ k},(x”')[j^_^ < oo. What is left is to show that 

m = mi{k]tk = T} is finite. To see this, without loss of generality assume that Ci > 1. 
Observe that 1/2 < + Vpiz)[t^_^^t^y < 2max(Hi(o)[4^_^,t^], for 

each k, which implies that (1/4)^ < + ^p(^)[4-i,tfc]> k gN. Consequently, 

^ m 

^ (^i(a)r + M^)t) < 00 , (3.9) 

k=l 

which completes the proof of (13.7p . 
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Theorem 3.3 Assume (F), (G) and that there exists h € such that I < h < u 

and Vp{h)T < oo, T € M"*". Let {{x"',k^)} denote the sequenee of Pieard’s iterations defined 
by iS. 5|) . Then for every T G R+, 

Vp{x^ — x)t —^ 0 and Vp{k^ — k)T 0, 

where (x, k) is a unique solution of 

Proof. Step 1. Convergence of Picard’s iteration. Fix T G M"*". Since X^ = Xq ^ = xo, 
applying Corollary 12.41 we get 

Vp(x^ - = Vp(x^ - x^-^fi 

<{d + l)Vp{[ f{s,x'fz'^) - f{s,x^z'^)das+ [ g{s, x'fZ^) - g{s, x'fz'^) dzs)t 
Jo Jo 

<{d + l)Vp{ [ fis, x:z^) - fis, 

Jo 

+ {d+ l)Vp{ [ g{s, x”r^) - g{s, x'fZ^) dzs)t 
Jo 

for t G [0,T]. By (|3.7I) . sup^<'r |x”| < for n € N, where N = sup„ Vp(x'^)T- Therefore 

Vpi [ fis,x^Z^) - f{s,x^z‘^)das)t < LNYfiafi sup\x^~^ - 
Jf) S<t 

< LNVi{a)tVp{x^-^ - x^-^)t 

and by ()3.2p . 

Vpi [ gis, x'fz'^) - gis, x^Z^) dzs)t < Cp^rVrigi-.,x'^~^) - g{-,x'^~‘^))tVpiz)f, 

Jo 

where r = (p/aAr)V(l//3). To estimate the right hand-side of the last inequality we will use 
the following lemma. 

Lemma 3.4 If g : M+ x M'’* — )> M satisfies (G) then for any x,y G B(M+,R'’*), T,N G R"*" 
such that Vpix)T < oo, Vpiy)T < oo, supj< 2 ’ \xt\, supj< 2 ’ \yt\ < N and r = (p/aN) V (1//3) we 


have 


Vrigiu x) - gi-,y))T < C^Yfix - y)T + Cn sup \xt - ?/t| (T^ + Vp{x)‘^^ + Vpiy) 

t<T ^ 

Proof of the lemma. For every t,s G [0,T], 

\g{t, xt) - git, yt) - gis, Xs) + gis, ?/s)| 


ojv 

T 


Vxgit, Oxt + (1 - 0)yt)ixt - yt) - "^xgis, 9xs + (1 - d)ys)ixs - ys)d6 


Hence 


\git, Xt) - git, yt) - gis, Xs) + gis, 2 /s)| 


: [ \Vxgit,6xt + il-9)yt)\\ixt-yt-Xs + ys)\d9 

Jo 

+ f \'Vxgit,9xt + il - 9)yt) -Vxgis,9xs + il - 9)ys)\\ixs - ys)\d9 


< C^lxt -yt-Xs + ys\ + Cn\xs - + \xt - + \yt - ?/s|“^). 
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Applying this estimate to each pair t = ti, s = ti-i from an arbitrary partition 0 = to < < 

■ ■ ■ <tn = T of [0, T] and using Minkowski’s inequality, we obtain the desired result. □ 
By Lemma 13.41 for i, j = 1,..., d we have 

< C^Vrix^-^-x^-\ 

+ Cn sup |xr 1 - xr^i + Vpix^-^)t'^ + Vpix^-^)t^), 

s<t ^ ^ 

which implies that 

d 

*d=i 

< {C^fWrix^-^ - x'^-% 

+ {CnY^ sviv\x'^-^ -x'^-^\U + Vp{x'^-^)'^^ +Vp{x'^-^)^ a. 

s<t ^ ^ 

From the above estimates, (I32D and fact that 

sup |x"-^ - x”-2| < l/p(x"-^ - x"-2)t, t € M+ 

S<t 

we conclude that there exists D > 0 depending only on Cp^r-, C^, Cn, un, f3, Ln, T and d 
such that for every n G N, 

Vpix'^ - x”-i)t < D{Viia)t + Vp(z)t)Vp(x^-^ - x^-% (3.10) 

Set ti = inf{t > 0; {D{Vi{a)t + Vp{z)t)) > ^} AT and observe that by induction, 

Vp{x^ - x”-i)t,_ < 2-(”-i)'Fp(x^ - x°)ti_, n € N. 

Thus {x”} is a Cauchy sequence in the space of cadlag functions on [0, ti) with the p-variation 
norm. Therefore there is a cadlag function x such that Vp{x'^ — x)t^_ —^ 0. This implies 
that Vp{jQ{f{s,x^_)- f{s,Xs-))das)ti- —^ 0 and Vp{f^{g{s,x^_) - g{s,Xs-)) dzs)ti- —^ 0, 
and hence that there exists a cadlag function k such that Vp{k^ — k)t^- —>■ 0 and (x, k) is a 
solution of (13.ip on the interval [0,ti). If we set 


Xt^ = max(min(xii_ + f{ti,Xti-)Aat^ + g{ti,Xti-)Azt^,Ut^),lti) 

and kti = kt^_ + Axt^ - {f{ti,xti-)Aat^ + g{ti,xt^-)Azt^) then by Remark [2T](b), (x,/c) is 
a solution of (|3.ip on the closed interval Moreover, 

A =max(min(x))^_ +/(ti, xJ\T^)Aati + g{ti,x^A)^^A^'^A),ki) 

—5- max(min(xti_ + f{ti,xt^-)Aat^ + g{ti,xt^-)Azt^,ut^),lt^) = xt^, 

which implies that i^(x”' — x)t^ —> 0 and Vp{kT — k)t^ —> 0. It is easy to see that we can 
apply the arguments used above to the interval [ti,t2] with t2 = inf{t > 0 ; T>(Vi(a)[t^ + 

^(^)[ti,t)) ^ and then to intervals [^2)^3]) [^3)^4])--- Since hiPa)^ < 00 and Vp{z)T < 

00, in finitely many steps we are able to construct the solution (x, k) of ( 13.111 on the whole 
interval [0, T] and to show that Vp{x"' — x)t — > 0 and Vp{kT — k)T —0. 
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Step 2. Uniqueness of solutions of (|3.1I1 . Assume that there exists two solutions 
and Let ti be defined as in Step 1. Using arguments from Step 1 we show that 

Vp{x^ — < \Vp{x^ — , which implies that x^ = on [0,ti). Since by Remark 

I2.ir bl we know that 

= max(min(x^^_ +/(ti,x^^_)Aaq + g{ti,xl^_)/^zt^,ut^),lt^), j = 1,2, 

it is clear that x^ = on the closed interval [0, ti]. Applying the above argument to intervals 
[ti,t 2 ], [^ 2 ,^ 3 ], ... we show in finitely many steps that x^ = x^ on [0,T] for every T G M+. □ 


4 Discrete-time approximation and stability of solutions 

We assume that l,h,u € are such that I < h < u, Vp{h)T < 00 for T G M"''. Let 

xo G [^0)^o] and a G D(R+,M), 2; G D(R+,M'^) be such that Ui(a)'r < 00 and Vp{z)T < 00, 
T G R+. 

Set Xq = Xo, /cq = 0 and 

^y1k+i)/n = fik/n, a:fc/„)(o(fc+i)/n - afc/n) + g{k/n, x"/^)(z(fc+i)/„ - z^/n), 

< x^k+i)/n =max(min(x^/^ + Ay”^^)/^,U(fc+i)/„),/(fc+i)/„), (4.1) 

t.n _ I ( n _ ™)i 'i _ A.,n 

'^(fc+l)/n — '^k ^ -^k/n) ^y(fc+l)/n 

and xf = x^^^, t G [k/n, {k + l)/n), A: G N U {0}. Since 

^'(k+l)/n = nC(fc+i)/„ {Xk/n + ^y(k+l)/n) ) A: G N U { 0 }, 

where ^C(k+i)/-a denotes the projection on the set = [l(k+i)/n,U(k+i)/n\, dE]) is the 

well known Euler scheme for dsn (see, e.g., @7]). It is also an analogue of the so-called 
“catching-up” algorithm introduced by Moreau to prove the existence of a solution of (II.Ih 
(see, e.g., [2]). 


Theorem 4.1 Let {(x”, A:”)} he a sequence of approximations defined by (EZP- If f, 9 satisfy 
(F), (G) and moreover f is continuous then 

{x^,k^,r,u^) ^ {x,k,l,u) inB(R+,R^'^), (4.2) 

where If = Ik/m uf = Uk/m t G [A:/n, (A; + l)/n), A: G N U {0} and {x,k) is a unique solution 

of 


Proof. FixT G R^ and set hr = Vi{a)T+ Vp{z)T+ Vp{h)T+ supf<^rpmax{\Alt\,\Aut\). First 
we show that if (x, k) satisfies (13.ip then 

yp{x)T < D where D depends only on d,xo, L,C^’'^, fi andbT- (4.3) 

By Corollary 12.51 for any t <T, 


Vpix)t < {d + l)Vpiy)t + dVp{h)t 


< (d + 1) 


\xo\ + yp{ f{s,Xs-)das)t + yp{ g{s,Xs-)dzs)t 
Jo Jo 


+ dypih)t. 
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We have VpiJ^ f{s, Xs-) das)t < Vi{a)t sup^<t |/(s,Xs_)| < LVi{a)t {l + Vp{x)t) and, by ([32]) 
and (13.41) . 




( / g{s,Xs-)dzs)t < Cp^p\/{i/p)Vpy(^i/i3){g{-,x))tVp{z)t 

Jo 

<Cp^p^ii/fS)C^'^{l + Vp{x)t)Vp{z)t. 


Hence there is Cq > 0 depending only on d, xq, Vp{h)T and Ci > 0 depending on d, L, /3, 
such that 

Vp{x)t < Co + Ci(l + Vp{x)t){Vi{a)t + Vp(z)t). 

Set ti = inf{t; Ci{Vi{a)t + Vp(z)t) > ^} AT. By the above, 

^p{x)[o,ti) < Co + - + 2^(^)[o,ti)) 
which implies that t^(x)[o^t^) < 2(Co + 1/2). Since by (13.3p . 


< \f{ti,xt^-)Aat^ \ + \g{ti,xt^-)Azt.^ \ + max(|A/tJ, |Ani^|) 

< Lr(l + \xti-\)\Aati\ + C/3,r(l + |xti_|)|AztJ + max(|AZiJ, |AtttJ), 

it is clear that t^(x)[o,ti] < D where D depends only on d, xq, C, C^C^ ^ and bx- If we 
set tk = inf{t > 4-i;Ci(Vi(a)[t^_^^t] + AT, k = 2,3,..., then for by 

the same arguments Vpix)it^_j^,tk] — ^ where D is depending only on d,xt^_.^,L,C^'"^, fd 
and bx- Since m = inf{A:;tfc = T} is bounded (similarly to (13.91) one can check that m < 
4 P(Vi(o)t + Vp{z)x)), this completes the proof of (14.3p . 

Now set a/ = a^/^, z^ = Zk/n, = hk/n, Pt = k/n, t £ [k/n, {k + l)/n), k gNU {0}. It 
is an elementary check that 

r {x^,k^) = ESP{y^,l,u), 

\ where = xq + /g /(p”_, x"_) da^ + g(p^_,x^_) dz^, t £ M+. 

Clearly, for any n € N, 


Hi(o")r < Hi(a)T, Vp{z^)x < V),(z)rand< Vp{h)x. (4.4) 

Combining (|4.3p with (|4.4p we get 

supt^(x"')'r < oo, T £ M^. (4.5) 

n 

Let xl"'^ = Xfc/„, k[^'^ = ky-j.^, t £ [k/n, {k + l)/n), k £ NU{0}, denote the discretization of the 
solution {x,k). By using [22l Chapter 3, Proposition 6.5] and [25l Chapter VI, Proposition 
2.2] one can check that {x^^'> ,a^, z"',P,u^, p^) —^ {x,a, z,l,u, I) in B(M+, R^'^+^), where 
Is = s, s £ M'*'. From this and an easy extension of [27\ Proposition 2.9] to functions with 
bounded p-variation it follows that 

[r = X 0 + ^ fip^-,x^;]) dal/ + ^ g{p/_,x^:}) dz/,P, u^) 

—>{y = XQ + J f{s,Xs-)das + J g{s,Xs-)dzs,l,v^ in B(M+,M^'^). 
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By the above and (12.211 


(x”, P, y”, r, •u’^) —s> (x, /c, y, I, u) in 




(4.6) 


where (x”,fe”) = ESP{y^^vT'), n € N. Moreover, analysis similar to that in the proof of 
(j4.3p shows that 


sup 1 ^(x”)t < oo, T € 


(4.7) 


By (j4.6|l and [25l Chapter VI, Proposition 2.2], (x"^,x^”^) —> (x,x) in ©(M'*',which 
implies that 

0, T G 


I -n M I 

sup |Xj — xl \ 
t<T 


Combining the above convergence with (j4.7|l and the fact that Vp{x^"'^)T < Vp{x)T < oo, for 
every e > 0 we obtain 

Vp+,{x^ - x("))t < Osc(x" - x(’^))y"^/^^+"^Vp(x" - ^0, T € M+, (4.8) 

where Osc(x)r = sup^ |xi — x^j. Fix e > 0. By Corollary 12.41 for any n G N and t <T, 


Vp+,{x^-x4 




<{d+l)Vp+,{j 

f{p4,x4)- 

f{p4,x^:^)d4+ f 

g{p4,x4)-gip4,x^4)dz^)t 

Jo 


Jo 


<{d+l)Vp+,i[ 

fip4,4-)- 

fip4^x^r-)d4 + [ 

g{p4,4-)- ff(p4 > ) dz^)t 

Jo 


Jo 


+ (d + l)Vpj^^{ 

[ f{p4,x4) 

- f{p4,4-)d4 + 

[ g{p4,x4)-9ip4,4-)dz^[ 

o 

'0 

J 

'o 

tTI,! I rn,2 

— • 





From (14.7p . the estimates from the proof of Theorem l3.3l and the fact that Ip+e(u)r < Vp{v)T 
one can deduce that there is I? > 0 such that for any n G N, 

4'" < DiV^ia^T + Vp{z^)T)VpU^^ - x (’^)) t . 

This together with (14.4p and (14.8p shows that lim„_>.oo 4^ = 0. The same arguments and 
(HSI) show that there is D > 0 such that for every t <T, 

4’^ < D{Viia'^)t + Vp{z4t)VpUx" - 

If we set ti = inf{t; D{Vi{a)t + Vp(z)t) > 1/2} A T then by (|4.4I) . 

Vp+,{x^ - x^4- < 4^ + ^VpUx^ - x4,. , 

which implies that Vp+e{x'^ — x"')t^_ ^ 0. This and (14.811 imply that Vp+e{x'^ — > 0. 

Note that 

x)} = max(min(x)}_ + /(ti,x)}_)Aa)} + y(ti, x)}_)Az[), u/J,//J, n G N. 
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If ti is a nonrational number then = AzJ^ = = Art”^ = Ax|”^ = 0. Hence 

and which implies that 

^0. (4.9) 

In case ti rational, set I = {n; there is k such that ti = k/n} and observe that if n € / then 
Aa”^ = Aaj^, = Azt^, Al^^ = A/^^, Au^^ = Aut^ and Axj”^ = Axt^. Consequently, 

lim„_^oo,ng/ = xt-^ and lim„_,.oo,ng/ ^ = xt^. On the other hand, if n ^ / then 

and = x|”l , which completes the proof of (I4.9h in case of rational ti. By the same 
arguments in finitely many steps we show that 

Therefore sup^^rp —>■ 0, from which we deduce that sup^^rp \ k‘^ — 0, which 

together with (|4.6p completes the proof of ()4.2I) . □ 

Corollary 4.2 Under the assumptions of Theorem \4-l[ for any T G M"*", 

^n<T ~ ^ ^ ^n<T ~ ^ 

where {x, k) is a unique solution of id. 1\) . 

Proof. By (14.2p and [25l Chapter VI, Proposition 2.2], 

^ {x,x,k,k) inD(R+,R^'^). 

This implies that x” — 0 in B(R''',R'^) and ^ 0 in B(R''',R'^), which is 

equivalent to ()4.10l) . □ 

Theorem 4.3 Assume (F), (G) and that there exists h G B(R+,R'^) such that I < h < u 

and Vp{h)T < oo for T G R"*". For e > 0 let Iq < x^ < uq and let fe,ge be functions satisfying 

(F), (G) with constants L,P,C^ not depending on e. If {x^,k^) denotes a solution of LI. 1\) 

with xo, f,g replaced by Xq, fe,gt and Xq —>■ xq, fe —>■ f, ge —>g then for every T G R^, 

/c /c 

I^(x'^ — x)t 0 and Vp{k'^ — k)T —>■ 0, 
where (x, k) is a unique solution of LI. 1\) . 


Proof. First observe that by (I4.3|l . 

supi^(x'^)'r < oo, T G 
£>0 

Fix T G R+. By Corollary 12.41 for every t G [0,T], 


(4.11) 


Vp(x" - x)t < (d + l)(^|xo - xol + Vp{ j fe{s,xl_) - f{s,Xs-)das)t 
9eis,xl_) - g{s,Xs-)dzs)t'^ 

< 


{d + l)(\x^o - Xol + Vp{J^ if, - f){s, xl_) das)t + Vp{J^ {g^ - g){s, x^_) dzs)t 
+ '^pij fis,xl_)-f{s,Xs-)das)t + Vp{J g{s,xl_) - g{s,Xs-)dzs)t'^ 

= (d + 1) (^|xo — Xol + . 
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Set N = sup^ Vp{x^)T, which is finite by to (|4.11l) Clearly, supj< 2 ’ \xl\ < N and 
= Vp{ [ {fe-f){s,xl_)das)T < sup \fe{t,v) - f{t,v)\Vi{a)T — 

Jo veB(0,N) 

t&[0,T] 

Since p < 2, there exists 7 € (1 — 1/p,/3 A (1/p)). Therefore by (j3.2p . (j3.4h and (I4.11h there 
is C > 0 depending only on 7 , L, /3, a,nd Vp{z)T such that 

4 ’^ = Vp{ [ {ge - g){s,xl_) dzs)T < Cp^i/^Vi/^{{ge - g){-,x^))TVp{z)T 
Jo 

< <4,i/7(Osc((pe -£/)(•, ((56 -p)(-,a;^))r^^^^^^^'^ 

+ \{9e- g){0,X^Q)\)Vp{z)T 

<C sup {\g^{t,v) - + \ge{0,v) - g{0,v)\). 

vGB{0,N) 
te[o,T] 

Consequently, 

+ (4.12) 

Using once again (14.lip and estimates from the proof of Theorem 13.31 we check that there is 
H > 0 such that for every t <T, 

{d + l){lf + if) < D{Vi{a)t + Vp{z)t)Vp{x^ - x)t. 

Similarly to the proof of Theorem 13.31 we set ti = inf{f; Zl(Ui(a)t + Vp{z)t) > 1/2} A T. 
Observe that 

Up(x" - x)ti- <id+ l)(|xo - xol + if + if) + ^Up(x" - x)ti- . 

From this and (j4.12p we deduce that Vp{x’^ — x)t^_ —> 0. Using arguments from the proof of 
Theorem 13.31 we show that this implies that Vp{x'^ — x)ti —^ 0. Applying this argument to 

(finitely many) intervals we prove the theorem. □ 

5 Applications to stochastic processes 

In this section we apply our deterministic results to SDEs with reflecting boundary condition. 
Let (Tt), P) be a hltered probability space and let A be an (J)) adapted process with 

trajectories in D(M+,M), Z,L,U,H be (J/) adapted processes with trajectories in B(R+,M'^) 
such that L < H < U and P{Vi{A)t < 00 ) = 1 , P{Vp{Z)T < 00 ) = 1 , P{Vp{H)T < 00 ) = 1 
for every T G M"*". Note that Z need not be a semimartingale. However, it is a Dirichlet 
process and a p-semimartigale in the sense considered in m and [29l [3^ . 

Definition 5.1 Let L <U and Xq be an Pq measurable random vector such that Lq < Xq < 
Uq. We say that a pair {X,K) of (T)) adapted processes with trajectories in ©(M+jM*^) such 
that P{Vp{X)T < 00 ) = 1 for r € R+ is a strong solution of (jl.ip if {X, K) = ESP{Y, L, U), 
where ^ ^ 

Yt = Xo+ [ fis,Xs-)dA,+ [ g{s,X,_)dZ„ f G R+. 

Jo Jo 


15 















Theorem 5.2 Assume (F), (G). Let L, U, H be {J't) adapted processes with trajectories in 
such that L < H <U and P{Vp{H)T < oo) = 1 . If Xq is an Xq measurable random 
vector such that Lq < Xq < Uq then lil-l]) has a unique strong solution {X,K). Moreover, 
if we define {(X"',iF”')} to be a sequence of Picard’s iterations for i.e. {X^,K^) = 

ESP{Xo,L,U) and for each n e N, = Xo + f{s, X^-^) dAs + g{s, X^-^) dZ^ and 
(X’^, iF”) = ESP{Y^, L, U) then for every T G R+, 

Vp{X^ - X)t ^ Q, P-a.s. and Vp{K^ - K)t ^ t), P-a.s. 


Proof. From Theorem 15.51 we deduce that for every cj G there exists a unique solution 
{X(ui) , K (lv)) = ESP{Y (w), L{uj), U{u!)) and for every T G R"'', 

Pp(X’"(w)-X( w))t ^ 0 , P-a.s., Pp(iF’"(a;) - iF(a;))T ^ 0 , P-a.s. 


Since for each n G N the pair (X”',iF"') is {Pt) adapted, the pair of limit processes {X,K) is 
{Ft) adapted as well, which completes the proof. □ 

Let be a fractional Brownian motion (fBm) with Hurst index H > 1/2 and let 
a : R"*" —>■ R is a measurable function such that 

Ikllfi/ff •= ( / < oo, T G R^. 

'^[O.T] Jo 

One can observe that the process a a dBf is a centered Gaussian process with 

continuous trajectories. Moreover, by |38l Theorem 1.1], for every r > 0, 

-<r < 

for all 0 < ti <t 2 - Hence for any subdivision vr = {0 = to < ■ ■ ■ < = F} of [0, T] we have 


Y.{E\Z^^ - < (C(l, Y. 

i=l i=l 



aa\^/^ds) = {C{l,H)f/^\\at{% . 


Therefore from [26l Theorem 3.2] it follows that if p > 1/Ft then P{Vp{Z^)T < oo) = 1 for 
T G R"*" (note also that Z^ is a Dirichlet process from the class studied in |14jL To 

approximate Z^ one can use the methods developed in [49] . 

We now show how to apply Theorem 15.21 and our approximation results of Section 4 
to fractional SDEs with reflecting boundary condition. Let B^ = {B^I ,... ,B^'‘^), where 
are independent fractional Brownian motions, and let Z^ = [Z^I ,..., Z^'‘^), 
where = f al dB^’^ with F : R"*" R such that |](T®|| i/h < oo for T > 0, z = 1,..., d. 

^[0,T] 

Let a : R"*" —>■ R be a continuous function with locally bounded variation. We consider 
fractional SDEs of the form 

Xt = Xty+ f f{s,Xa)daa+ f g{s,Xa)dZ^ + Kt, fGR+. (5.1) 

Jo Jo 

Glearly, (15.ip generalizes classical fractional SDEs driven by B^ and is a particular case of 

m- 
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For any re G N we set 


X^ = Xli^, K^ = Kli^, tG[A:/re,(A; + l)/re), A:GNU{0}, (5.2) 

where Xq = Xq, Kq = 0 and 

' ^^il+wn = /(*/«. - at/n) + 3{k/n, - Z«J. 

^(kj^l)/n — ^k/n^ y^{k+l)/n ^k/n) ^^(k+l)/n’ 

Corollary 5.3 Assume (F), (G). Let L, U, H be (J^) adapted processes with continuous 
trajectories such that L < H < U and P{Vp{H)T < oo) = 1 for T G If Xq is an Xq 
measurable random vector such that Lq < Xq < Uq then 15.il) has a unique strong solution 
{X,K). Moreover, i/{(X”, iF"')} is a sequence of approximation defined by 15.Dl) and 15.,'?)) 
then for every T G M"*", 

sup \Xff -Xt\^0, P-a.s., sup \Kf - Kt\ 0, P-a.s. (5.4) 

t<T t<T 

Proof. It suffices to apply Theorem 15.21 and Theorem 14.11 The uniform convergence follows 
from the fact that if a is a continuous function and Z^,L,U have continuous trajectories 
then also the solution (X, K) has continuous trajectories. □ 

Note that in the case where L, U may have jumps Theorem 14.II implies weaker then (|5.4I) 
convergence. Namely, we then have 

{X^,K^) iX,K), P-a.s. in D(M+,m2"'). 

6 Proof of Theorem 12.21 

We follow the proof of [23l Theorem 2.1]. 

Step 1. We assume additionally that y^,y‘^ and the barriers l,u are step functions of the 
form 

yl = Y/, k = Li, ut = Ui, t€[ti_i,ti), i = l,...,n-l 

and yl = Yi, It = L^, ut = Un, t G [tn-i,tn = T], j = 1,2, for some partition 0 = 
io < < ■ ■ ■ < tn = T of the interval [0,T]. By Remark I2.11 bh ki = Ki,te 

i = l,...,re - 1, k{ = Ki, t G [tn-i,tn = T], j = 1,2, where = 0 and = 

max(min(iF/_p Pj - Yi),Li - Y^), i = 1,... ,n, j = 1,2. Clearly 

L^-Yi <Kl <Ui-Yi, i = l,...,n, 3 = 1,2. (6.1) 

Without loss of generality we may and will assume that Vp{k\ — kg)T > 0. Hence there 
exists i such that iF/ ^ or Kf ^ Later on, without loss of generality we will 

assume that for any i = 1 ,..., re — 1, 

Kl^Kl, or Kf^Kl, (6.2) 

(If (j6.2l) does not hold then we set vq = 0, 

Vk = inf{i > Vk-i-,K} 7^ iF/.^oriFf / iFf_ J A re, A: = 1,. .. ,re. 
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re = inf{A:;tifc = re}, yj = h = L^,,, Ut = C4j^, t € [tvk_,ty^) for A; = 1,...,re - 1, 
y{ = Y^, It = Ln, Ut = Un, for t G = T], j = 1,2. Then (|0]) holds true for the 

functions y^, = ESP{y^ J,u), j = 1,2, and moreover, Vp{k^ — k‘^)T = Vpik^ — k‘^)T 

and Vp{y^ - y‘^)T < Vp{y^ -y‘^)T- Consequently, 

Vp{k^ - P)r < Vp{f - fh Vp{k^ - k^)T < Vp{y^ - y^h-) 


It is clear that there exist numbers Q = io < ii < ... < im = n such that 

m 

- k^)T = Y. K^’i - YJ - (6.3) 

k=l 


and 

{Kl -KU-{kI-KI_J^ 0, k = l,..., m. (6.4) 

Hence, if m > 2 then for k = 2,... ,m we have 


UK-, - Y-U - U<L, - - Y-,)- K, -KK < «. 


Indeed, if (16.5p is not satisfied then by (|6.4p . 


(6.5) 


1(^-3 - - (K-t - + l(K - KJ - - Kjr 

< i(K-t - K-J - (Ki-t - Klj + (^ - Ki-t) -iK- ^-jr 

= \{Kl-KU-{Kl-Kl_J\P, 

which contradicts ([ 631 ). 

We will show that there exists 0 = io < rf < ii (resp. 0 = io < < ii) such that if 

- Kf^ < 0 (resp. - Kf^ > 0) then 


and for k = 2, 


- y3 


<Kl-Kl 


(resp.i^i -Kl <y3-b;V) 


m there exist < rujvY < such that 


r_A — r„A 




( 6 . 6 ) 


(6.7) 


Fix A: = I,..., rre. Set = maxji < : Kl = Ut — Y^ oiKl = Li — Y?}, j = 1, 2, with the 

convention that max0 = 0. By (16.21) . or r^. = Without loss of generality we may 

and will assume that = ik- Then we have three cases: 


(a) 

Y - 

II 

{L^,- 

-y.M 

ik> 


-b;|) (or 

k 

Y- 

II 


y-M - 

iu^l - 

-rp) 

(b) 

Y- 

*4 = 

^ik — 

Yl^ and rl = 

0 (or Kl 

-Y 


- Yl and rl 

= 0), 


(c) 

Y- 


{Lik- 

-'•'L 

-Kl 

-y|) (or 

^k 

Y- 

Kl = 

iUr,- 

r.;)- 

(L,2 - 

-1|)) 

By p6.ip in 

all the 

cases 











Kl- 

- k'^ = 

ik 

' -bjfe - 

-d- 

■fits 

Lik - Yl 

-Kl 

+ Kl 

-{Lik 

- ft) 

= Y^ - 

Ik 

- Y^ 

*fe’ 
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which implies that we can put = ik- In order to find we consider the cases (a), (b), (c) 
separately. 

In case (a), if = ifc then 

A'i - Kj = i.. - r.; - (i.. - Y,l ) = Y,l- F,; 

and we put = ik- H ’>'1 < ik then we set r* = max{i < i^ : Kj = Ui — Y^} M r‘f.. Observe 
that = ■ ■ ■ = Since for r* < v < ik, = max(iirt,_i, Ly — Yy), it follows 

by (16.21) that 

Kl. - <---<Kl-Kl. (6.8) 

From this it also follows that F),* — 1),* < — Kf^. Indeed, if r* > (resp. r* = r'^ ) then 

Ky^, = Ur* — Yy„ (resp. = Lr* — F)* ) and by ()6.1I) . 

Kl - Kl > Kl. - > Ur* - f;1 - Ur* + Yyl = - Yy\, 

(resp. — Ky* > Ly* — F),* — Ly* + = ^r* ~ Ir*)- What is left is to put 

= r* and show that ik-i < r*. This is obvious if A: = 1, so assume that k > 2 and 
ik-i > r*. By ()6.8p . — {Kf^ — > 0. From this ()6.5p it follows that 

~ ^ik- 2 ^ ~ ~ ^ik- 2 ^ ^ which together with p6.8p implies that ik -2 < r*. 

Using once again p6.8p we see that {Ky* — — {Ky* — < 0. Consequently, 

0>(AL.-A-y)-(AL.-Ay) 

> (A.L, - UJ - (kU - Klj + (At - A.L,) - (K% - Kl_^) 

= {K}., - - {K% - 

and 

\{Kl_^ - KlJ - {Kl_^ - Kl_^)\P < \{Kl. - KlJ - {Kl - kI_1\U (6.9) 

Similarly, 

^<iKl-Kl_l-{Kl-Kl_l 

< {Kl - Kl_l - {Kl - Kl_l - {Kl - Kl_l + {Kl - kI_ 1 
= {Kl - Kl) - {Kl - Kl), 

which implies that 

\{Kl - Kl_l) - {Kl - Kl_l)\P <\{Kl- Kl) - {Kl - kI)\U (6.10) 

Combining p6.9p with p6.10p we obtain 

- KlJ - {Kl_^ - KlJl + \{Kl - Kl_l -{Kl- Kl_l)\^ 

< \{Kl - Kl_l - {Kl - Kl_l\P + |(iF^ - Kl) - {Kl - kI)\p, 

which contradicts (16.3p and completes the proof of the fact that ik-i < r*. Consequently, in 
case (a) we put = r*. 
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In case (b) (resp. (c)) we set r* = max{z < ik : = Ui — Y^} (resp. r* = max{z < ik : 

Kf = Li — oiKI = Ui — Y-^}). For r* < u < 4 we have = max(iF^_^, — Y^) and 

— Y^). As in case (a) we conclude from this and ()6.2p that that 

By the argument used in case (a) we also show that ik-i < r*. Moreover, if r* > 0 and 
= Lr* — Y^* (resp. A^* = Ur* — then by (|6.1I) . 

Kl -Kl> Kl. - {Lr* - Y,l) > {Lr* - Y,\) - {Lr* - y4) = >;* - 

(resp. Kl - Kl > {Ur* - Y^) - K^ > {Ur* - Y^) - {Ur* - Y^.) = Y^r - y/*)- Therefore 

we put = r*. Since /c = 1 if r* = 0, the proof of (16.61) and (16.7p is complete. 

Now observe that by (16.71) . if Kl — Kl > Kl__^ — A?^for some k = 2,... ,m then 


0 < - Kl) - {Kl_^ - KlJ < {Yl - Yl) - {Yl_^ - YlJ, 

which implies that 

\iKl -Kl)- (At, - Kl_X < I (tv - tv) - (t^, - t^,)r• 

Similarly, if A/^ - Kl < A^, - 

\iKl - Kl) - (At, - KlJ\^ < Kt - t) - (t^, - tv_,)r- 

In case A: = 1, if K} — Kf > 0 then 

L\ L\ 

\{Kl-Kl)- {Kl - Kl)\P = \Kl -Klr< Itv - tvT 

and if Kf — Kf < 0 then 
61 61 

\{Kl - Kl) - {Kl - At r < It - tfl"- 

Putting together (|6.1ip ~ ()6.14D we conclude that 


( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

(6.14) 


E ml - A) - (A-. - Urn" < Y, - ill I”+E iwi -Yj- (uL. - uL.)r. 


k=l 


k =2 


where fk = t ^ < ik, k = 1,... ,m. Hence 

Vp{k^ - klr = Vp{k^ - klr 


i/p 


= (Ei(A-A)-(A-.-'A.: 

k=l 

m ^ I 

< (iiV - Yir + E 1(4 - 4) - (4-. - YLif) 

k=2 

m ^ I 

< |yo - 2/01 + (I (2/t - yl ,) - (t_, -1,_,) I'’) 


k=l 
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for some partition 0 = tfQ < < ■ ■ ■ < < T, which proves the theorem in the case of 

step functions and step barriers l,u. 

Step 2. The general case. 

Let {^”'1 and {n”} be sequences of discretizations of I and u, 

respectively, i.e. y]''^ = ylj^, y^’” = yl/^,lt = 4/n> < = u^/n t € [k/n,{k + l)/n), k € 
NU{0}. By [25l Chapter VI, Proposition 2.2], {y^’^,y‘^’^,l^,u^) —)• {y^,y‘^,l,u) inD(M+,M^). 
Let = ESP{yk^J^^u^), n e N, j = 1,2. By {kk^,k‘^’^,yk^,y^’^) 

{k^, k'^, y^, y'^) in B(R+,M^), which implies that 

kl,n _ J^2n - ykl _f.2 D(M+, M). (6.15) 

By Step 1, for n G N and T G M+ we have Vp{k^’'^ — k‘^’'^)T < i^(y^’”' — y^’”)T- Clearly, 
Vp{y^’^ — y‘^'^)T < yp{y^ — y'^)T, n G N, T G M^. From this and (j6.15p it follows that for 
every T G R'*' such that Ak^ = Ak"^ = Ay^ = Ay|, = 0, 

Vp{k^ - e)T < liminf - k^’^)T < sup Vp(yi’" - y^’-)^ < Vp{y^ - y%. 

n—>-oo ^ n 


To obtain the desired result for arbitrary T G R^ we use right continuity of Vp{k 
Vp{y^ - 2 /^)- 


1 


k‘^) and 
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